Black holes in the brane world : Time symmetric initial data 
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We numerically construct time-symmetric initial data sets of a black hole in the Randall-Sundrum 
brane world model, assuming that the black hole is spherical on the brane. We find that the apparent 
horizon is cigar-shaped in the 5D spacetime. 



I. INTRODUCTION 

Motivated by Horava-Witten model the so called 
brane world model has been actively investigated 
Among several models, a simple, but very attractive 
model was recently proposed by Randall and Sundrum 
|^,||. According to their scenario, we are living in a 4D 
domain wall in 5D bulk spacetime. The noteworthy fea- 
tures of their model are that in the linearized theory, the 
conventional gravity can be recovered on the brane |^-|^] 
and that a homogeneous, isotropic universe can be sim- 
ply described if we consider a 4D domain wall moving in 
the 5D Schwarzschild-anti de Sitter spacetime ||^ . 

One of the most non-linear objects in the theory of 
gravity is a black hole, which should be also investi- 
gated to understand the nature of the models in strong 
fields. However, because of the complexity of the equa- 
tions, any realistic, exact solutions for black holes have 
not been discovered in the brane world model, even with 
help of numerical computation so far. We only know 
that the effective 4D gravitational equation on the brane 
is different from the Einstein equation Q (see Appendix 
A), so that the static solution for a non-rotating black 
hole should not be identical with the 4D Schwarzschild 
solution. Indeed, a linear perturbation analysis [^] 
shows that a solution of gravitational field outside self- 
gravitating bodies on the brane is slightly different from 
the 4D Schwarzschild solution. Chamblin et al. con- 
jecture that the topology of black hole event horizons 
would be spherical with the cigar-shaped surface in the 
5D spacetime. However, nothing has been clarified sub- 
stantially. 

In this paper, as a first step toward self-consistent stud- 
ies for black holes in the brane world, we numerically 
compute a black hole space using a time symmetric initial 
value formulation; namely we solve the 5D Einstein equa- 
tion only on a spacelike 4D hypersurface. Thus, the black 
hole obtained here is not static nor the exact solution for 
the 5D Einstein equation, implying that we cannot iden- 
tify the event horizon. However, we can investigate the 
property of the horizon determining the apparent horizon 



which could give us an insight on the black hole in the 
brane world. We focus on the Randall-Sundrum's second 
model Q , and assume that the black hole is spherical on 
the brane, but the shape of the horizon is non-trivial in 
the bulk. We will determine the apparent horizon on the 
brane and show that the black hole is cigar-shaped as 
conjectured in pO|. 



II. FORMULATION AND RESULTS 

We consider time symmetric, spacelike hypersurfaces, 
St, in the brane world model assuming the vanishing 
extrinsic curvature; i.e., 



((5"-l-i^r)(4)V„i, = 0, 



(2.1) 



where t^ is the unit normal timelike vector to Ej and 
'^"'^Vq is the covariant derivative with respect to the 4D 
metric on St. In this case, the momentum constraint is 
satisfied trivially, and the equation of the Hamiltonian 
constraint becomes 



16^G5(A + (5)T^^tA.i-)^ 



(2.2) 



where ^^^i? is the Ricci scalar on St, and G'5(= k^/Stt), 
A and '■^•'T^jy denote the gravitational constant, negative 
cosmological constant, and energy- momentum tensor in 
5D spacetime [cf., Eq. (Al)]. We choose the line element 
on St in the form 



dl^ 



1 r 



i^dz"^ + ^P^{dr^ -^r^dn) 



(2.3) 



where £ = -y/— K5A/6, z (> 1) denotes the coordinate 
orthogonal to the brane and r (> 0) is the radial coordi- 
nate on the brane. We assume that the brane is located 
at z = 1. Note that we simply choose this line element 
for convenience of the analysis. In this paper, we fo- 
cus on a black hole which is spherical on the brane, i.e., 
ip = il}{r, z). Then, the explicit form of the Hamiltonian 
constraint in the bulk (for z > 1) is written in the form 



1 



2,,^ 3 



K-2 



(2.4) 



where ' = 9/9r, and ^^V^^ is the energy-momentum ten- 
sor in the bulk, which is introduced for numerical conve- 
nience. 

Equation (2.4) is an elliptic type equation and should 
be solved imposing boundary conditions at z = 1, z 3> 1, 
r = 0, and r ^ £. The boundary condition at z = 1 
is derived from Israel's junction condition as (see 
Appendix A for the derivation) 



(2.5) 



The boundary conditions at z ^ 1 and r ^ £ are ob- 
tained from the linear perturbation analysis (see Ap- 
pendix B). For r ^ £ and r > £z, it becomes 



MGa 



2r 



(2.6) 



where G4 = G5/£, M is the gravitational mass on the 
brane, and R = (2/3) ^/^f. For z > 1, 



1+4- 



3G4M 



Rz 



-3/2 



(2.7) 



To determine the existence of a black hole, we search 
for the apparent horizon. Here, we determine two hori- 
zons One is defined to be the spherical two-surface 
on the brane on which the expansion of the null geodesic 
congruence confined on the brane is zero M], i.e., 



93 = |j(2V'' + i^)=0. 



(2.8) 



The other is the apparent horizon in full 4D space, which 
is defined with respect to the n ull g eodesic congruence in 
full 5D spacetime and satisfies 



= 0, 



(2.9) 



where is a unit normal to the surface of the apparent 
horizon. Explicit equation for determining this apparent 
horizon is shown in Appendix C. 

The procedure of numerical analysis is as follows. 
First, we artificially put the matter of ph = ^^V^ixi^i" > 
in the bulk. This method is employed because we do 
not have to consider the inner boundary condition of 
black holes with this treatment. As long as ph is confined 
around the brane and inside the horizon, it does not sig- 
nificantly affe ct th e geometry outside the horizon. Then, 
we solve Eq. (^), and try to find the apparent horizon 
both on the brane and in the bulk. When the distribu- 
tion of Ph is sufficiently compact, the apparent horizons 
exist. It should be noted that two horizons do not coin- 
cidently appear. In some cases, the apparent horizon on 
the brane exists although that in the bulk does not. 
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FIG. 1. Location of the apparent horizons on the brane 
(filled circle) and in the 4D space (solid line) . Artificial matter 
is confined in the region shown by the dashed line. 
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FIG. 2. Profile of -i/; — 1 on the brane (solid line). Location 
of the apparent horizon on the brane is shown. The dashed 
line denotes — 1 = M /2r 

Here, we show one example of numerical results. We 
set G4 = 1. In this example, an artificial matter is put for 
< r < Q.2R and 1 < z < 1.2. Equation is solved 
using a uniform grid with grid size 1200 x 1200 for r and 
z directions, which covers a domain with < r / R < 17.1 
and 1 < z < 18.1. In this case, the gravitational mass 
on the brane is M ~ 0.29i?, and both apparent horizons 
on the brane and in the bulk exist. We note that the re- 
sults are essentially the same for 0.25 < M/R < 0.5. In 
Fig. 1, we show the location of apparent horizons in the 
bulk and on the brane. The apparent horizon in the bulk 
is apparently cigar-shaped. Due to this cigar-shape the 
circumferential radius of the apparent horizon is different 
depending on the choice of the circumference in the bulk. 
In Fig. 2, we show that the profile of — 1 on the brane. 
For r ^ R, t/j — 1 behaves as M/2r, implying that the 
solution approximately agrees with that in the 4D Ein- 
stein gravity, i.e., the bulk effect is small. However, the 
existence of the bulk is significant for r ~ i? as expected. 
Indeed, ip — 1 deviates from M/2r with decreasing r. This 
effect is in particular important for the location and area 
of the apparent horizon on the brane: In the case of 4D 
gravity without bulk, the apparent horizon is located at 
7" AH — M/2 with the area Aah — IGttM^. However, in 
the brane world model, they take different values in gen- 
eral. (In this example, tah — 0.9M and Aar ^ 88. 6M^, 
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and the coefficients converge to well-know 4D values (0.5 
and IGtt) with increasing M , implying that the effect of 
the existence of the bulk becomes less important.) 



III. SUMMARY 

We numerically computed time symmetric initial data 
sets of a black hole in the brane world model, assuming 
that the black hole is spherical on the brane. As has 
been expected, the black hole (apparent horizon) is cigar- 
shaped in the bulk |l0| . 

We remind that we only present time symmetric ini- 
tial data of a black hole space. This implies that the 
black hole is not static and will evolve to other state 
with time evolution. The quantitative features of the 
final fate could be different from the present result. Self- 
consistent analysis for static black holes should be carried 
out for future to obtain a definite answer with regard to 
black holes in the brane world. However, we believe that 
the present result provides us a guideline for such future 
works. 
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cases in which ^^^T^j^ = 0. Using (4-1-1) formalism, the 
effective 4D equation on the brane has the form 



(A3) 



where ^'^-'G^jy is the 4D Einstein tensor on the brane, 
A4^^4(^A+^4X^) and = ^^^C^p.^n^n^, (A4) 
where '"^'^C^pua is 5D Weyl tensor. In the above, for sim- 



plicity, we set '■^V^i, = 0. Equation (A3) implies that 



we can consider E'^t^ as the effective source term of the 
4D Einstein equation on the brane, and as long as E^^^ 
is not vanishing, the geometry on the brane is different 
from that in the 4D gravity even in the vacuum case. 
Only for very special case such as for the black string 
solution (10|P|], i^p^ = holds. 

From Eq. (|A3|) , we find that the Minkowski spacetime 
is realized on the brane when E^^, — and A4 = 0. In 
this paper, we set A4 = to focus on asymptotically flat 
brane. Then, the junction condition at x = is rewritten 



to K, 



In the case when we choose the line 



elem ent as Eq. (2^), the junction condition reduces to 
Eq. (O). 
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APPENDIX A: THE ESSENCE OF THE BRANE 
WORLD 

We briefly review the covariant formalism of the brane 
world 1^. For the matter source of the 5D Einstein 
equation, ^^^G^^ = kK^^^T^^ - A'^'g^,,), we choose the 
energy-momentum tensor as 



p,u\ 



' pv J 



(Al) 



where % = £lnz, A is the tension of the brane, is the 
induced metric on the brane, and ^^^T^ji/ is the energy 
momentum tensor on the brane. Due to the singular 
source at x = and the Z2 symmetry, we can derive the 
Israel's junction condition at x = as 



j^pv — g'^s'^yM"^ 2 v 3 



(A2) 



where K^j^^, = q'^q^D^nx, and and are the covari- 
ant derivative with respect to qpi, and the unit spacelike 
normal vector to the brane. In the text, we consider the 



APPENDIX B: ASYMPTOTIC BOUNDARY 
CONDITIONS 

To specify the boundary condition at infinities, we in- 
vestigate the linearized equation of Eq. (2.4): 



(Bl) 



where ip — I + ip and (p I. We can obtain the formal 
solution with aid of the Green function G{x, z; a;', z') as 



2T:Gd / d^x' dz' G{x, z]x' , z')ph{x' , z'). (B2) 



Assuming that ph is non-zero only in the small region 
around the brane, we can derive the relevant Green func- 
tion as ||] 



G(x, z; x' , z') — — 



(27r)3 
1 

JO 

Go + Gkk, 



'tk-{x— 5< 



dm 



k2 



(B3) 



where Um{z) is the mode function given from the Bessel 
functions J„ and Nn as 



Ura{z) = Z' 



mR? 
2t 



Ji{mR)N2{mRz) - Ni{mR)J2{mRz) 
X — ^ -, (B4) 



V(Ji(mi?))2 + (7Vi(mi?))2 



3 



where R = {2/3y^'^£. Gq and Gkk are the Green func- 
tion of zero and KK modes, respectively. From Eq. (B2) 
we can derive the asymptotic boundary conditions shown 
in the text. 



APPENDIX C: APPARENT HORIZON IN THE 
BULK 

We derive the equation for the apparent horizon in the 
bulk. After we perform the coordinate transformation 
from (r, z) to {x, 9) as z = 1 + a;| cos9\ and r = ^a;sin6', 
the surface of the apparent horizon is denoted by a: = 
h{9). Then, the non-zero components of Si is written as 



C and S0 



-Ch, 



(CI) 



At 9 = 7r/2, the boundary condition is imposed as h,g = 
0. 

Note that in the limit 9 ^ tt/2 (i.e., on the brane), Eq. 
is written in the form 



(C6) 



where we use hg = and the relation dgip = D = dxC = 
0. Note that the equation which the apparent horizon on 
the brane satisfies is Ad^-^j/ip + 2//i = [cf., Eq. (2.8)]. 
Thus, unless Cp'h/d9'^ = h a,t 9 = tt/2, the apparent hori- 
zon on the brane cannot coincide with that in 4D space. 
Note that the black string solution exceptionally 
satisfies (Fh/d9^ ^ h at 9 ^ n/2. 



where C[= ip'^C/{l + x\cos9\)] is a normalization con- 
stant calculated from s*Si = 1, and h^g = dh/d9. Then, 
the equation for h can be written to the following ordi- 
nary differential equation of second order 



d^ 



, dx'ib 3 dxC 

^JCY_ _| ^ _x_ 

Tp h{l + h\cos9\) C 



X sm 



V'^ cos^6l- (1 - V )sin6'cos6'-^ 



+h- 



■0 



h sin 9 



l + h\ cos9\ 



+ 2cot9 + D 



X 1^(1 - V' ) sin 6* cos 9 - (cos^ 9 + ?// sin"" 9) 

+4:ip^ dxipicos'^ 9 + h^^ sm9cos9h,0) 

+h^^(l - tp"^) sin 9 cos 9h^e 

+h-\l - 7/-*) cos(26i) - ih-^ sin 9 cos ^V'^^eV' 



where 
D 



+{(l-V/)sin(20) -4sin2 0^^^dg^}-^ 



-C^[(l - V/){1 - /i-^/ifJ sine cos 6* 
-/i"^(l - V/)cos(26')/i,e 
+2i:^dgij{cos 9 + h~^ sin 9h^gf] 



(C2) 



(C3) 



Eq. ( |C2| ) is solved imposing boundary conditions at = 
and 7r/2. In the limit — > 0, we impose the following 
boundary condition. 



h = ho + /i20^ 



where /12 is evaluated at x 
following equation; 



8dx^ 



0(9^), 
- ho and 

dxC 



(C4) 
from the 



Tp hoil + ho) C 



ho 



(C5) 
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